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ABSTRACT 

An  inverse  Finite  Element  Method  (iFEM)  is  presented  for  beam  and  frame 
structures.  The  method  is  aimed  at  reconstructing  the  complete  displacement  field 
starting  from  in  situ  measurements  of  surface  strains.  A  laboratory  experiment  is 
conducted  on  a  simple  cantilever  beam  subjected  to  various  static  loadings. 
Experimentally  measured  strains  are  used  within  a  single-element  iFEM  model  to 
assess  the  efficiency  and  predictive  capability  of  the  approach  with  respect  to 
uncertainties  and  measurement  errors  that  unavoidably  affect  real  structures. 

INTRODUCTION 

Real-time  reconstruction  of  the  deformed  shape  of  a  structure  is  a  key  technology 
for  Structural  Health  Monitoring  (SHM)  [1].  The  computation  of  the  displacement 
field  is  perfonned  on  the  basis  of  in  situ  strain  data  computed  in  real  time  by  a  network 
of  strain  sensors.  This  inverse  problem  is  commonly  referred  to  as  shape  sensing. 

Few  approaches  in  the  literature  lead  to  practical  algorithms  applicable  for  real¬ 
time  shape  sensing  of  aerospace  structures  (for  a  more  complete  review  refer  to  [1,2]). 
Ko  [3]  proposed  a  one-dimensional  scheme  based  on  classical  beam  theory,  in  which 
the  deflection  and  cross-sectional  twist  angle  of  an  aircraft  wing  are  determined  by 
establishing  the  position  of  the  wing’s  neutral  axis,  using  a  representative  finite 
element  solution  of  a  wing.  The  axial  surface-strain  measurements  are  specified  at  the 
nodes  of  a  piecewise  approximation  of  the  idealized  beam;  the  method  is  assessed 
using  the  strains  from  a  FEM  solution.  Maint^on  and  co-workers,  [4,5],  developed  a 
finite  element  formulation,  seeking  the  solution  for  displacements  and  loads 
simultaneously.  Starting  from  measured  displacements  and  strains  and  requiring  a 
priori  knowledge  of  the  material  properties  and  a  subset  of  applied  loading,  the 
formulation  results  in  the  number  of  unknowns  that  is  three  times  the  number  of 
degrees-of-freedom  (dof  s)  in  the  finite  element  discretization.  In  [5]  the  effect  of 
modeling  errors  is  investigated.  The  approach  was  also  assessed  on  a  damaged  beam 
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using  experimentally  measured  data.  Tessler  and  Spangler  [1]  developed  an  inverse 
Finite  Element  Method  (iFEM)  for  shear-deformable  plate  and  shell  structures.  The 
formulation  is  based  on  a  least-squares  variational  principle  that  allows  the  full-field 
reconstruction  of  the  three-dimensional  displacement  vector  from  measured  surface 
strains.  An  experimental  assessment  of  the  method  was  presented  in  [6,7],  where  the 
deformed  shape  of  a  slender  beam  is  reconstructed  by  an  iFEM  shell  model  using 
Fiber-Bragg  Grating  (FBG)  strain  measurements.  Recently,  Gherlone  [8,9]  presented 
an  iFEM  formulation  for  shear-deformable  beam  and  frame  structures.  Beam  and 
frame  shape-sensing  analyses  were  performed  for  static  [10]  and  dynamic  [2]  loadings 
using  strain  data  from  high-fidelity  FEM  models. 

The  paper  presents  a  brief  review  of  the  iFEM  variational  formulation  that 
incorporates  experimentally  measured  strains  within  a  simple  inverse  beam-frame 
element.  The  element  is  based  on  Timoshenko  beam  theory  which  includes  the  axial, 
bending,  torsional  and  transverse  shear  deformations.  A  laboratory  experiment  is  then 
discussed  in  which  a  thin-walled  cantilevered  beam  is  tested  using  a  set  of  static 
loadings.  The  experimentally  measured  strain-rosette  data  are  then  used  within  a 
single  inverse-beam  element  to  model  the  entire  beam.  The  modeling  accuracy  of  the 
iFEM-reconstructed  deformed  shape  of  the  beam  is  assessed  by  comparing  the 
predicted  displacements  with  those  measured  experimentally  by  the  displacement 
transducers  distributed  along  the  beam’s  span. 

INVERSE  BEAM-FRAME  FINITE  ELEMENT  FORMULATION 

Consider  an  isotropic,  straight  beam-frame  structural  member  of  Young’s  modulus 
E ,  shear  modulus  G  ,  and  Poisson  ratio  v  (Figure  1(a)).  The  structural  member  is 
referred  to  a  Cartesian  coordinate  system  (x,  y,z),  where  x  is  positioned  along  the 

centroidal  and  shear  axis,  and  y  and  z  are  the  cross-section’s  principal  inertial  axes. 
The  frame  member  has  length  L,  cross-sectional  area  A,  area  moments  of  inertia  with 
respect  to  the  y  -  and  z  -axis,  Iy  and  I_ ,  and  polar  moment  of  inertia  Ip  =  I  +  /. . 

Consistent  with  the  hypotheses  of  Timoshenko  beam  theory  (each  cross-section 
remains  flat  and  rigid  with  respect  to  thickness-stretch  deformations  along  the  y-  and 
z-axis)  [11]  and  neglecting  axial  warping  due  to  torsion,  the  displacement  field  is 

u.x  {x,  y,  z)  =  u(x)  +  z0y  (x)  -  y07  (x) 

uy(x,y,z)  =  v(x)~  z6x  ( x) ;  uz(x,y,z)  =  w(x)  +  yOx  ( x) 

where  ux,  u  ,  and  a  are  the  displacements  along  the  x -,  y  -,  and  z -axes 
respectively;  u  ,  v ,  and  w  are  the  corresponding  average  displacements;  0X ,  <9V ,  and 
6  are  the  rotations  about  the  three  coordinate  axes.  The  kinematic  variables , 

u  =  jt,  v,  w,  <9  ,  0y,  6  j  and  their  positive  orientations  are  shown  in  Figure 
1(a).  The  displacement  field  (1)  gives  rise  to  the  linear  strains 


ex  (x,  y,  z)  =  e1  (x)  +  z  e2  (x)  +  y  e3  (x) 

yxz (x,  y)  =  eA{x)  +  y e6(x);  yxy(x,z)  =  e5 (x)  ~ze6 (x) 


(2) 


-IT 

e6  \  are  given  as 


where  the  strain  measures  e(u)  =  [en  e2,  e3,  e4,  e5, 


e\  ~  U,x  ’  e2  ~  @y,x  ’  g3  —  @z,x  ’  eA  ~  W,x  ’  g5  —  V ,x  ^z  ’  —  ^x,x  (^) 


The  iFEM  formulation  reconstructs  the  deformed  structural  shape  by  minimizing  a 
weighted  least-squares  functional  <£>  containing  the  strain  measures  obtained  by  in- 
situ  strain  sensors,  e£ ,  and  e(u)  defined  by  Eqs.  (3),  i.e., 

<b(u)  =  ||e(u)-ef||  (4) 

The  ‘Flu)  functional  is  then  discretized  by  ^-continuous  finite  elements  with  the 
displacements  u  defined  by 


u(x)-u,!(x)  =  N(x)u''  (5) 

where  N(x)  denotes  ^-continuous  shape  functions  and  u  the  nodal  dofs. 
Consequently,  the  total  least-squares  functional  is  a  sum  of  the  N  individual  element 
contributions,  <F'’(u/'),  i.e.,  <F  =  ^ A  <J>‘’ .  Accounting  for  the  axial  stretching, 
bending,  twisting,  and  transverse  shearing,  the  element  functional  is  given  by  [2] 

<I>V)  =  XA^  (6) 

k= 1 


where  I  =  [/l”,  A"  Ae,  A3Ae,  A",  A°,  A"  4'  ] ,  with  A]'  (k  =  denoting 

dimensionless  weight  coefficients,  and  Ae  the  element’s  cross-sectional  area.  The  six 
components  of  the  element  functional  are  given  as  the  Euclidean  norms 

2 

4>; --£[<•,  W-ef]  (*  =  1 . 6)  (7) 

11  1=1 

where  Le  is  the  element  length,  n  is  the  number  of  strain  sensors  within  an  element, 
and  xt  ( 0  <  xt  <  If)  are  the  positions  at  which  the  strain  sensors  are  located,  with  the 
superscript  a  denoting  the  strain  measures  that  are  computed  from  the  strain  sensor 
values  (experimental  values)  at  the  xf  locations.  Invoking  Eqs.  (2)  and  (5),  the  analytic 
element-level  strain  measures  are  expressed  in  matrix  form  as 

e(u)  =  B(x)u£’  (8) 

where  the  matrix  B(x)  contains  the  derivatives  of  the  shape  functions  N(x) . 
Substituting  Eq.  (8)  into  Eq.  (7)  and  minimizing  the  element  functional  with  respect  to 
uc  results  in  the  element  matrix  equation  keue  =ke,  where  ke  depends  only  on  the 


measurement  locations  xt ,  whereas  the  vector  f e  depends  on  the  experimentally 

measured  strains  [9].  Assembly  of  the  element  contributions  for  the  total  problem 
discretization,  while  taking  into  account  transformations  from  the  element  (local)  to 
the  global  coordinate  system,  and  upon  specifying  the  problem  boundary  conditions, 
results  in  a  non  singular  system  of  matrix  equations  of  the  form  KU  =  F .  The  solution 
of  these  equations  for  U  is  efficient:  the  K  matrix  is  inverted  only  once,  since  it  is 
independent  of  the  measured  strains  and  remains  unchanged  for  a  given  distribution  of 
strain  sensors.  The  F  vector,  however,  is  dependent  on  the  measured  strain  values; 
thus,  at  any  strain  measurement  update  during  deformation,  the  matrix-vector 
multiplication  provides  the  solution  for  the  unknown  nodal  displacement  dofs 
contained  in  U. 


(a)  (b) 

Figure  1.  (a)  Beam  geometry  and  kinematic  variables,  (b)  inverse  finite  element  geometry  and  dof  s. 

For  frame  structural  elements  loaded  only  by  forces  and  moments  at  the  end 
points,  the  strain  measures  can  be  shown  to  exhibit  the  following  span-wise 
distributions:  ex,  eA ,  e5,  and  e6  are  constant,  whereas  e2  and  e3  are  linear  [9].  From 

Eq.  (3),  it  is  deduced  that  u  and  6X  are  linear,  9  and  0  parabolic,  v  and  w  cubic. 
Thus,  the  following  interpolations  are  adopted  [2,9,10]  (also  refer  to  Figure  1(b)): 


«(#)=LC(#K  «,(#)=!  *!"(#)«„ 

1=1,2  /=1,2 

V(£)  =  E  (#)v,  -  I  Nf<  (f )  A,;  e,  (f)=  E  A  (#)«*  (9) 

*=1,2  7=1, r, 2  7=1, r, 2 

»■(#)= E^(#)w,+  E  AT™ (#)»„:  *,(#)=  E 

1=1,2  7=1,  r,  2  7=1, r,  2 

where  E,  =  (2x  /  Zf  - 1)  e  [-1,  l]  is  a  non-dimensional  axial  coordinate;  the  subscripts  1, 
r  and  2  denote  positions  along  the  beam  length  at  £  =  —  1,  <^  =  0,  and  £  =  +l, 
respectively;  L'1’  (^)  (z  =  1, 2)  are  linear  Lagrange  polynomials;  Lf'  (g)  (7  =  1,  r, 2) 
are  quadratic  Lagrange  polynomials.  The  cubic  N(p  (g)  (7  =  1, r, 2)  polynomials  are 
obtained  from  standard  cubic  Lagrange  polynomials  by  enforcing  the  transverse  shear 
strain  measures  (e4  and  e5 )  to  be  constant  along  the  element  (refer  to  [2]  for  the 


expression  of  N |3)  (g) .)  The  element  has  fourteen  dof  s:  six  at  each  end  node  plus 
the  rotations  0ym  and  0,m  at  the  mid- span.  The  two  internal  rotation  dofs  are 
condensed  out  statically  to  achieve  a  two-node,  twelve  dof  element. 


v<9 

(b) 


x 

-> 


(a) 


Figure  2.  (a)  Orthogonal  and  cylindrical  coordinate  systems,  (b)  strain  gauge  location  and  coordinates. 

A  key  step  in  the  formulation  is  to  compute  the  strain  measures,  e£,  from  the 
experimentally  measured  surface  strains.  We  shall  restrict  the  present  analysis  to  the 
beam-frame  members  with  circular  cross-sections  only,  and  employ  the  cylindrical 
coordinate  system  (0,x,  r)  shown  in  Figure  2(a).  A  strain  gauge  is  placed  on  the 

external  surface  ( r  =  Rext ),  at  x  =  xi ,  and  is  oriented  along  the  6  and  [3  angles  (refer 
to  Figure  2(b)).  The  relation  between  the  measured  strain  e2  and  the  six  strain 
measures  at  x  =  x;  is  [8,9] 


+e4  {xt)cpspcg-e5  (x)c^s>y0+e6  (-t)vA* 


where  cg  =  cos  6 ,  sg  =  sin  9 ,  cp  =  cos  jB  ,  and  sp  =  sin  jB  .  For  the  case  of  end-node 
forces  and  moments,  the  constant  distributions  of  ex,  e4,  <?5 ,  and  eb  and  linear  e2  and 
e3  can  be  determined  requiring  a  total  of  eight  strain  measurements.  This  number  may 
be  reduced  to  six  by  invoking  the  equilibrium  equations  of  Timoshenko  beam  theory 
which  relate  the  bending  moments  ( Af  ,  M  _ )  to  the  transverse  shear  forces  ( Qy ,  Q  ) 


where  k2y  and  k:  denote  the  shear  correction  factors  [9,11]. 


EXPERIMENTAL  RESULTS 


The  application  of  iFEM  for  beam  and  frame  structures  has  been  carried  out  using 
a  laboratory  experiment  conducted  at  the  AERMEC  Laboratory  of  the  Aeronautics  and 
Space  Engineering  Department  of  Politecnico  di  Torino.  The  test  article,  shown  in 


Figure  3(a),  is  a  thin-walled  circular  cross-section  beam  having  thickness  s= 2  mm, 
external  radius  Rexi= 40  mm  and  length  £= 800  mm.  The  material  is  a  6060  aluminum 
alloy  (£’=61,922  MPa,  v=0.33).  The  beam  has  been  mounted  on  a  test-bed  in  a 
cantilevered  configuration  with  one  end  clamped  between  two  couples  of  iron  blocks 
locked  together  by  twelve  bolted  connections.  Four  different  load  cases  have  been 
considered  (Figure  3(b)):  (i)  tip  vertical  force,  (ii)  tip  horizontal  force,  (iii)  tip  force 
inclined  at  30°  with  respect  to  the  horizontal  y-axis  of  the  cross-section,  and  (iv)  tip 
vertical  force  applied  at  a  distance  6=3 00  mm  from  the  center  of  the  cross-section.  The 
loading  has  been  achieved  by  placing  several  weights  on  a  cradle  (for  a  total  weight  of 
£=26.83  Kg).  For  cases  (i),  (ii),  and  (iii),  the  cradle  was  linked  to  a  screw  at  the  center 
of  the  beam  tip  cross-section  while;  for  load  case  (iv),  the  cradle  was  suspended  to  a 
proper  lever  arm  provided  by  a  thick  plate  embedded  at  the  beam  tip  (Figure  3(a)). 


(a)  (b) 

Figure  3.  (a)  Experimental  set-up  (A:  clamping  system,  B:  lever  arm  for  load  case  (iv),  C:  loading  system, 
D:  LVDT  transducers,  E:  aluminum  beam),  (b)  load  cases. 


Nine  stacked  strain  rosettes  have  been  placed  at  three  different  stations  along  the  beam 
(x=L/ 3,  LI 2,  2L/3);  for  each  station,  the  three  rosettes  are  placed  at  6=- 120°,  0°,  +120°, 
respectively.  Each  rosette  has  three  strain  gauge  measuring  at  /?= 0°,  45°,  90°  (Figure 
2(b)).  Table  I  shows  the  strain  gauge  configurations  used.  In  the  notation,  the  first 
number  indicates  the  number  of  strain  measurements,  and  the  letter  “E”  signifies  the 
use  of  Eqs.  (11).  LVDT  measurements  have  also  been  taken  at  different  locations 
along  the  beam  length  (Figure  3(a))  and  used  to  assess  the  iFEM-recovered  deflections 
and  rotations.  For  all  load  cases,  a  single  inverse  element  is  used  to  model  the  entire 
beam,  using  =1  (k  -  1, ...,  6) ,  see  Eq.  (6). 

Table  II  summarizes  the  experimental  and  iFEM  results.  For  each  loading,  the 
measured  deflections  and  rotations  are  shown,  together  with  the  percent  difference 
with  the  corresponding  iFEM  predictions. 


TABLE  I.  STRAIN  GAUGE  DISTRIBUTIONS 
(ORIENTATIONS  ARE  EXPRESSED  IN  DEGREES). 


Notation 

Orientation  (6,  ft)  at  x=U 3 

Orientation  (6,  J3)  at  x=LI2 

Orientation  (  0.f3)  at  x=2LI3 

6Ea 

- 

(-120,0),  (-120,45),  (0,0), 
(0,45),  (120,0),  (120,45) 

- 

6Eb 

(-120.0),  (0,0),  (120,0) 

- 

(-120,0),  (0,45),  (120,0) 

8a 

(-120,0),  (-120,45),  (0,0), 
(120,45) 

- 

(0,0),  (0,45),  (120,0), 
(120,45) 

8b 

(-120,0) 

(-120,0),  (-120,45),  (0,0), 
(0,45),  (120,0),  (120,45) 

(120,0) 

TABLE  II.  EXPERIMENTAL  AND  IFEM  RESULTS:  DEFLECTIONS  AND  ROTATIONS 
MEASURED  BY  LVDT  AND  PERCENT  DIFFERENCE  (IN  BOLD)  IN  IFEM  PREDICTIONS 
USING  A  SINGLE  INVERSE-ELEMENT  MODEL. 


Load 

case 

(i) 

Measured  displacements 
(mm) 

w(L/2) 

w(3U4) 

w(7U8 ) 

w(L) 

-0.648 

-1.300 

-1.635 

-1.985 

iFEM 

%  difference 

6Ea 

-9.48 

-7.82 

-5.11 

-3.09 

6Eb 

-2.85 

-2.45 

-0.36 

0.88 

8a 

-1.07 

-0.86 

0.01 

2.22 

8b 

-1.52 

-1.12 

0.92 

2.12 

Load 

case 

(ii) 

Measured  displacements 
(mm) 

v(U2) 

v(3U4) 

v(7U8) 

v(L) 

0.636 

1.272 

1.648 

1.994 

iFEM 

%  difference 

6Ea 

-4.73 

-3.70 

-4.34 

-2.61 

6Eb 

-4.00 

-3.59 

-4.58 

-3.24 

8a 

-5.24 

-3.98 

-4.62 

-2.95 

8b 

-2.87 

-1.99 

-2.87 

-1.42 

Load 

case 

(iii) 

Measured  displacements 
(mm) 

v(3U4) 

v(7U8 ) 

v(L) 

w(3U4) 

w(7U8) 

w(L) 

1.098 

1.414 

1.667 

0.644 

0.802 

0.963 

iFEM 

%  difference 

6Ea 

-8.81 

-8.44 

-3.85 

1.61 

4.66 

6.98 

6Eb 

-4.83 

-5.28 

-1.48 

-0.70 

2.40 

4.81 

8a 

-5.97 

-6.00 

-1.79 

-1.65 

1.68 

4.14 

8b 

-4.35 

-4.55 

-0.46 

1.93 

4.83 

6.75 

Load 

case 

(iv) 

Measured  displacements 
(mm)  and  rotations  (rad) 

w(L/2) 

w(3U4) 

w(7L/8) 

w(L) 

0X(L) 

-0.664 

-1.319 

-1.647 

-1.958 

3.5xl0"J 

iFEM 

%  difference 

6Ea 

-11.15 

-8.70 

-5.39 

-1.37 

-9.36 

6Eb 

-7.15 

-5.69 

-2.90 

0.50 

-5.23 

8a 

-4.10 

-3.12 

-0.58 

2.51 

-7.84 

8b 

-5.37 

-3.83 

-0.95 

2.55 

-9.36 

In  all  loading  cases  and  using  any  of  the  considered  strain-gauge  configurations,  the 
iFEM-predicted  tip  deflections  differ  from  the  measured  values  by  less  than  7%.  The 
deflections  computed  along  the  beam  span  exhibited  slightly  higher  differences, 
however,  not  exceeding  12%.  The  iFEM  model  that  used  eight  strain  gauges  generally 
produced  more  accurate  predictions  compared  to  the  six-strain  gauge  configuration, 
with  the  exception  of  case  (iv)  corresponding  to  6Eb  which  produced  a  slightly 
superior  prediction  for  the  tip  twist  rotation. 

The  results  of  this  first  experimental  assessment  of  iFEM  for  beams  confirm 
the  accuracy  and  efficiency  of  the  present  formulation,  particularly  taking  into  account 
that  only  a  single  inverse  element  was  used  to  model  the  beam. 


CONCLUSIONS 


An  inverse  finite  element  method  (iFEM)  has  been  presented  for  the  full-field 
reconstruction  of  the  Timoshenko-type  beam-frame  kinematics,  using  experimentally 
measured  strains  as  input  quantities  in  the  formulation.  The  approach  has  been 
demonstrated  on  a  cantilevered  beam  tested  under  several  static  loading  cases  in  a 
mechanics  laboratory.  Strains  measured  by  different  sets  of  strain  gauges,  attached  to 
the  beam’s  external  surface,  have  been  employed  in  the  inverse  analyses.  Using  only  a 
single-element  model,  the  iFEM-reconstructed  displacements  and  rotations  compared 
favorably  with  those  measured  experimentally.  The  effectiveness  of  the  iFEM 
modeling  has  been  shown  to  depend  on  the  layout  and  number  of  strain  gauges.  The 
results  of  the  present  effort  also  point  towards  the  possibility  of  determining  optimally 
distributed  locations  of  the  strain  gauges  in  order  to  achieve  further  improvements  in 
the  shape-sensing  predictions. 
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